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ABSTRACT 

We present an extended grid of mean three-dimensional (3D) spectra for low-mass, pure-hydrogen atmosphere 
DA white dwarfs (WDs). We use C05BOLD radiation-hydrodynamics 3D simulations covering T^ff = 6000- 
11,500 K and logg = 5-6.5 (g in cm s“^) to derive analytical functions to convert spectroscopically determined 
ID temperatures and surface gravities to 3D atmospheric parameters. Along with the previously published 3D 
models, the ID to 3D corrections are now available for essentially all known convective DA WDs (i.e., logg 
= 5-9). Eor low-mass WDs, the correction in temperature is relatively small (a few per cent at the most), but the 
surface gravities measured from the 3D models are lower by as much as 0.35 dex. We revisit the spectroscopic 
analysis of the extremely low-mass (ELM) WDs, and demonstrate that the 3D models largely resolve the 
discrepancies seen in the radius and mass measurements for relatively cool ELM WDs in eclipsing double WD 
and WD -i- milli-second pulsar binary systems. We also use the 3D corrections to revise the boundaries of the 
ZZ Ceti instability strip, including the recently found ELM pulsators. 

Keywords: convection - hydrodynamics - white dwarfs 

1. INTRODUCTION 

Recent three dimensional (3D) radiation hydrodynamical 
simulations of the Sun and main-sequence stars have led 
to significant improvements in our knowledge of stellar at- 
mospheres and spectroscopic abundance measurernents (e.g . 
lAsplund et alj|20()^ ICaffau et al.ll201 ll : IScott et ^l2015alm . 

The major differences between the ID and 3D model atmo¬ 
spheres result from the insuf ficient description of convection 
by the mixing length theory (lBohm-VitenselfT958l MET) that 
is used in the ID models. 

Convection becomes important in hydrogen-atmosphere 
(DA) white dwarfs (WDs) below about 15,000 K. Spec¬ 
troscopic analyzes of large samples of WDs with ID 
model atmospheres sh o w the so-called “high l o gg problem” 

(iBergeron et all 119901: lEisenstein et alJ [200d iKepler et alJ 
120071 iGianninas et all 1201 ll iKleinman et al.l 120131) , where 
the spectroscopically determined masses of WDs cooler than 
13,000 K are systematically higher than their hotter counter¬ 
parts by as much as 20%. There is no known evolution ef¬ 
fect or observational bias that should lead to such a change 
in mass for the cooler WDs, and a larger average mass is not 
observed for cool WDs with parallax measu rements. In a se¬ 
ries of papers. ITremblav et al.l (1201 lLl2013allbll3) presented the 
first 3D model atmosphere calculations for D A WDs using 
BOLD radiation-hydrodynamics code (lErevtag et alJ 
The predicted Balmer line profiles are significantly 
different in the 3D models, and the “high logg problem” es¬ 
sentially disappears in the spectroscopic analyzes based on 
the 3D models. The previously published model grid covers 
the effective temperature range of 6000-15,000 K and the sur¬ 
face gravity range of logg = 7-9. To study the effects of 3D 
models on the physical parameters of the extremely low-mass 
(ELM) WDs with logg < 7, here we extend our 3D model 
atmosphere grid down to logg = 5. 


the C O^ 
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ELM WDs are helium-core remnants and almost always 
found in short period binary systems. The ELM Survey 

(iBrown et al .1120101 : iKilic et al]l2010al IGianninas et al.ll2()14aL 

and references therein) has discovered binaries with period s 
as short as 12-20 min (iBrown et alJl201 ll iKilic et al.ll2014|) . 
Improving the constraints on the ELM WD physical parame¬ 
ters, including the spectroscopic mass and gravitational wave 
strain, is important for future gravitational wave missions 
in the milli-Hertz frequency range. In addition, all of the 
known eclipsing and/or tidally distorted double WD systems 
involve low-mass WDs. Studying the m ass-r adius relation 
for the se systems, IGianninas et all (1201 4a^ and iKaplan et alJ 
(1201 4ft find that the ID spectroscopic analysis underestimates 
the radius (or overestimates the surface gravity) for ELM 
WDs cooler than about 10,000 K. NETT 11748 is an excel¬ 
lent example, where the eclipses constrain the radius and sur- 
face gravity to R = 0.0423-0.0433 R© and logg = 6.32-6.38 
(IKaplan et al.ll2014ft, while the ID spectrosc opic analysis pre¬ 
dicts logg = 6.83 (IGianninas et al.ll2014aft . Therefore, the 
“high logg problem” seems to impact all convective DA WDs, 
including ELM WDs. 

We outline the details of our 3D model calculations in Sec¬ 
tion 2, and present the ID to 3D atmospheric parameter cor¬ 
rection functions for the extended range of logg = 5-9 in 
Section 3. We discuss the astrophysical implications, includ¬ 
ing the revised parameters for the ELM Survey sample, radius 
constraints from the eclipsing and tidally distorted systems as 
well as ELM WD -i- milli-second pulsar binaries, and the ZZ 
Ceti instability strip in Section 4. We conclude in Section 5. 

2. 3D MODEL ATMOSPHERES 

We have carried out 30 simulations of pure-hydrogen 
ELM WD atmospheres with CO^BOLD (lErevtag et alJ201^ . 
These 3D simulations in the range 5.0 < log? < 6.5 pro- 
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Figure 1. Surface gravity (logarithmic value) and mean Teff for our 
CO^BOLD 3D model atmospheres. The simulations were computed with 
a bottom boundary layer that is open (open circles) or closed (filled circles) 
to convective flows. 

vide a d irect extension to the grid presented in iTremblav et al.l 
(12013d) for 7.0 < logg < 9.0. The calculations also comple¬ 
ment t he CIFIST grid of CO^BOLD simulations for stars and 
giants (iLudwig et al.l l2009bt iCaffau et al.l 1201 Ih with 1.0 < 
logg < 4.5. Figure [1] illustrates the new ELM WD 3D sim¬ 
ulations in a HR-type diagram. Properties of the individual 
models, such as Teff (derived from the temporally and spa¬ 
tially averaged emergent stellar flux), logg, box dimensions, 
and computation time, can be found in Table[T] 

The numerical setup fo r the ELM WD simulat ions is the 
same as that described in iTremblav et al.l (I2013cl) for a grid 
of 150 X 150 X 150 point s. The prope r ties o f the boundary 
conditions are reported in lErevtag et al.l (120121 see Sect. 3.2). 
The warmest simulations, represented by filled symbols in 
Eigure [T] have a convective zone that is thinner than the ver¬ 
tical dimension of the atmosphere. In those cases, we use 
a bottom layer that is closed to convective flows with im¬ 
posed zero vertical velocities. In cooler models, the convec¬ 
tive flux is transported in and out of the domain through an 
open boundary. In all cases the lateral boundaries are peri¬ 
odic and the top boundary is open to material flows and radi- 


tables as those described in ITremblav et al. 

(2013d). 

[n par- 

ticular, we emnlov the iHummer & Mihalas 

(HU EOS, the 

Stark broadening profiles of Tremblav & 

lereeroni ( 

2009), 


and the quasi-molecular line opacity of I Allard et all ( 2004 ). 


The wavelength-depende nt opacities are s orted i n the 8 bin 
configuration discussed in iTremblav et al.l (l2013cl) . 


2.1. Physical Properties 

The ELM WD 3D simulations show physical properties that 
are very similar to those of the 3D models of higher mass 
WDs presented in earlier studies. ITremblav et al.l (12013bl) 
found that 3D effects on the thermal structure are well cor¬ 
related with the density at Rosseland optical depth tr = 1. 
Since the atmospheric density can be kept roughly constant 
by decreasing both logg and Teff, the ELM simulations have 
properties similar to those of C/O-core WDs but with a shift 
towards smaller Teff values. 

Figure |2] compares our simulation at ~7500 K and 
logg = 5.0 with the simulation at ~9500 K and logg = 7.0 


Table 1 

Grid of CO^BOLD 3D Model Atmosplieres for ELM WDs 


7’eff 

logg 

logz 

logx 

log'fR,min 

log tr, max 

log! 

loS ^adv 

(K) 

(cm s“^) 

(cm) 

(cm) 



(s) 

(s) 

5977 

5.00 

7.56 

8.38 

-6.58 

2.96 

4.00 

1.17 

6958 

5.00 

7.71 

8.46 

-6.08 

2.99 

4.00 

1.03 

7513 

5.00 

7.62 

8.33 

-5.64 

2.99 

4.00 

0.92 

8501 

5.00 

8.36 

8.64 

-7.77 

3.06 

3.30 

0.94 

9001 

5.00 

8.38 

8.64 

-5.46 

3.00 

3.30 

1.19 

9499 

5.00 

8.49 

8.68 

-8.29 

3.01 

3.30 

1.54 

6044 

5.50 

7.04 

7.66 

-6.98 

2.97 

3.50 

0.75 

6915 

5.50 

7.20 

7.73 

-6.58 

3.01 

3.50 

0.62 

7967 

5.50 

7.37 

7.98 

-6.57 

3.02 

3.50 

0.47 

8996 

5.50 

7.93 

8.18 

-7.47 

3.02 

2.90 

0.46 

9499 

5.50 

7.93 

8.18 

-6.12 

3.00 

2.90 

0.67 

9998 

5.50 

8.00 

8.22 

-6.75 

3.01 

2.90 

0.91 

10500 

5.50 

8.06 

8.29 

-7.23 

3.01 

2.90 

1.40 

5958 

6.00 

6.66 

7.20 

-7.21 

2.99 

3.00 

0.38 

6971 

6.00 

6.78 

7.28 

-6.99 

3.02 

3.00 

0.21 

7962 

6.00 

6.89 

7.34 

-6.82 

2.99 

3.00 

0.09 

8963 

6.00 

7.00 

7.42 

-5.76 

2.96 

3.00 

-0.08 

9491 

6.00 

7.32 

7.65 

-5.07 

3.00 

2.40 

-0.03 

9999 

6.00 

7.46 

7.71 

-7.34 

3.00 

2.40 

0.21 

10499 

6.00 

7.52 

7.75 

-5.75 

3.01 

2.40 

0.40 

11001 

6.00 

7.49 

7.75 

-5.00 

2.96 

2.40 

0.77 

5986 

6.50 

6.10 

6.57 

-7.01 

2.99 

2.50 

-0.03 

6922 

6.50 

6.20 

6.65 

-7.48 

3.00 

2.50 

-0.19 

7990 

6.50 

6.35 

6.84 

-6.76 

2.99 

2.50 

-0.31 

8907 

6.50 

6.51 

7.03 

—6.66 

3.00 

2.50 

-0.42 

9412 

6.50 

6.54 

6.95 

-5.29 

3.02 

2.50 

-0.49 

9989 

6.50 

6.95 

7.24 

-6.75 

3.02 

2.50 

-0.49 

10500 

6.50 

6.99 

7.24 

-5.80 

3.01 

2.05 

-0.32 

11000 

6.50 

7.05 

7.28 

-5.77 

3.01 

2.05 

-0.13 

11499 

6.50 

7.09 

7.33 

-5.83 

3.01 

2.05 

0.12 

6011“ 

5.00 

8.04 

8.08 

-7.93 

8.38 

5.00 

1.08 


“Stellar simulation witli Nuc/Nh = 8.51% and [M/H] = —4.0. 

Notes. All quantities were averaged over 12 snapsliots. r^ff is derived from 
the temporal and spatial average of the emergent flux, x (same as v) and z cor¬ 
respond to the geometrical dimensions of the box, and tr is the Rosseland 
optical depth averaged over constant geometrical depth, t is the computa¬ 
tion time in stellar' time and is the advective or turnover timescale at the 
geometrical depth that corresponds to (TR)x,y = 1. 


from ITremblav et aH (I2013cl) . Both simulations have a sim¬ 
ilar relative intensity contrast and granulation is visually 
comparable. Table |2] demonstrates that ELM WDs have a 
slightly larger maximum relative intensity contrast compared 
to higher mass degenerates. ELM WDs are very similar to 
main-sequence simula tions in that respect (see Figure 10 of 
ITremblav et al.ll^l3bl) . Furthermore, Table |2] describes the 
characteristic size of granulation and decay time. These quan¬ 
tities f ollow the same trends as those observed for higher mass 
WDs (ITremblav et al.ll2013bl) . in particular the relation be¬ 
tween the characteristic size of granulation to pressure scale 
height ratio and photospheric Mach number. 

Our low gravity WD simulations reach photospheric plasma 
conditions and characteristic geometrical scales that are sim¬ 
ilar to main-sequence stars, although the radii are quite dif¬ 
ferent. In order to better understand this link, we have com¬ 
puted a stellar model at Teff ~ 6000 K, logg = 5.0, = 
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Table 2 

Granulation Properties of ELM WDs 


Teff 

logg 

log Char. Size 

log//p 

lo§^decay 

<5/rms/{/> 

log Mach 

(K) 

(cm s“^) 

(cm) 

(cm) 

(S) 

(%) 


5977 

5.00 

7.50 

6.70 

1.53 

16.91 

-0.42 

6958 

5.00 

7.54 

6.74 

1.62 

23.95 

-0.26 

7513 

5.00 

7.57 

6.77 

1.65 

23.17 

-0.12 

8501 

5.00 

7.83 

6.86 

1.51 

18.34 

-0.07 

9001 

5.00 

7.69 

6.92 

1.45 

11.17 

-0.26 

9499 

5.00 

7.77 

6.99 

1.49 

3.80 

-0.58 

6044 

5.50 

6.78 

6.21 

1.12 

13.00 

-0.50 

6915 

5.50 

6.91 

6.25 

1.15 

21.22 

-0.34 

7967 

5.50 

7.08 

6.31 

1.17 

23.14 

-0.14 

8996 

5.50 

7.42 

6.41 

1.19 

18.83 

-0.06 

9499 

5.50 

7.27 

6.47 

0.98 

13.01 

-0.22 

9998 

5.50 

7.32 

6.53 

0.79 

6.93 

-0.42 

10500 

5.50 

7.46 

6.62 

0.80 

1.84 

-0.86 

5958 

6.00 

6.23 

5.71 

0.86 

8.02 

-0.62 

6971 

6.00 

6.40 

5.77 

0.80 

17.55 

-0.43 

7962 

6.00 

6.53 

5.82 

0.80 

22.37 

-0.27 

8963 

6.00 

6.63 

5.88 

0.80 

21.21 

-0.04 

9491 

6.00 

6.81 

5.94 

0.80 

19.55 

-0.05 

9999 

6.00 

6.82 

6.01 

0.58 

13.76 

-0.24 

10499 

6.00 

6.83 

6.07 

0.37 

7.65 

-0.39 

11001 

6.00 

6.81 

6.15 

0.64 

2.13 

-0.73 

5986 

6.50 

5.73 

5.21 

0.43 

5.13 

-0.71 

6922 

6.50 

5.89 

5.28 

0.37 

12.49 

-0.52 

7990 

6.50 

6.00 

5.33 

0.33 

19.63 

-0.36 

8907 

6.50 

6.13 

5.38 

0.37 

21.18 

-0.21 

9412 

6.50 

6.31 

5.42 

0.38 

20.45 

-0.10 

9989 

6.50 

6.47 

5.48 

0.30 

19.67 

-0.06 

10500 

6.50 

6.40 

5.55 

0.17 

16.19 

-0.19 

11000 

6.50 

6.35 

5.60 

0.00 

9.80 

-0.35 

11499 

6.50 

6.36 

5.68 

-0.20 

4.90 

-0.56 

6011“ 

5.00 

7.34 

6.60 

1.76 

17.79 

-0.38 


“Stellar simulation with Nhc/Nh = 8.51% and [M/H] = —4.0. 

Notes. All quantities were averaged over 250 snapshots except T^ff (see Ta¬ 
ble 1). The characteristic horizontal size of the granulation, decay time, and 
relative intensity contrast (5/rms/(f)) were c omputed from emergen t inten¬ 
sity snapshots using the definitions given in ITremblav et all j2013bl) . Both 
the pressure scale height (Hp) and the Mach number are computed at the 
geometrical depth that comesponds to {Tp)x,y = 1. 


8.51%, and a metallicity of [M/HI = -4.0 ba s ed on solar ele¬ 
mental abundance s from | G revesse & Sauvall (119981) with up¬ 
dated CNO values (lAsplundl2005h . The model is an extension 
of the CIEIST grid for main-sequence stars where logg < 4.5. 
The simulation is representative of the extremely metal-poor 
E-dwarfs with the currently highest derived surface gravities 
(iSnite et al.ll2012h . The properties of the simulation are de¬ 
scribed in Tables 1 and 2, while Eigure [3 shows the compar¬ 
ison with the pure-hydrogen ELM WD at Teff ^ 6000 K and 
logg = 5.0. The relative intensity contrast and size of the gran¬ 
ules are similar, although the stellar atmosphere is denser by 
58% at tr = 1 compared to the more opaque pure-hydrogen 
WD atmosphere. It implies that their properties are only ex¬ 
pected to be qualitatively similar. Nevertheless, it demon¬ 
strates that ELM WDs will provide an important connection 
between WDs and main-sequence stars in future studies of the 
3D effects on the determination of the atmospheric parameters 


T,„ = 7513 K, log g = 5.0, pure-H, 8 I„„ 3 /<I> = 23.17% 



0.5 1.0 1.5 2.0 
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5 10 15 20 25 
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Figure 2. Emergent bolometric intensity at the top of the computational box 
for the pure-hydrogen WD simulations at Feff = 7513 K, logg = 5.0 (top), 
and 9521 K, logg = 7.0 (bottom). The root-mean-square intensity contrast 
with respect to the mean intensity (5/ims/(7)) is identified above the panels. 
The length of the bar in the top right is ten times the pressure scale height at 
= 1 - 

from Baimer line profiles. 

We have computed mean 3D structures, hereafter (3D) 
structures, which are spatial and temporal averages of and 
P performed over surfaces of constant tr and for 12 random 
snapshots. We have employed those snapshots to define the 
Teff values of our simulations (Table 1). In Eigure |4l we 
compare the temperature structures of (3D) simulations and 
ID models at Teff ^ 7000 K in the range 5.0 < logg < 9.0. 
Convection is still fairly adiabatic in the photosphere for 
log^ > 7.0 in this Teff range, resulting in small differences 
between the (3D) and ID cases except for the optically thin 
layers where 3D overshoot causes a significant dynamic cool¬ 
ing. Density decreases for lower surface gravities and con¬ 
vection becomes less efficient in the photosphere. Eigure |4] 
demonstrates that differences between (3D) and ID structures 
become larger for ELM WDs at Teff 7000 K. The observed 
pattern is very si milar to the one a t consta nt logg and variable 
Teff in Eigure 7 of iTremblav et al.l ([20133) • 

3. 3D ATMOSPHERIC PARAMETERS 
3.1. Mean 3D Spectra 

We have computed spectra from the 3D simulations in or¬ 
der to study the 3D effects on the determination of atmo¬ 
spheric parameters. Eirst of all, we rely on mean spectra 
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= 6011 K, log g = 5.0, [M/H] = -4, bljd) = 17.79% 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

X (Mm) 


Figure 3. Similar to Figure[2|but for the emergent bolometric intensity of the 
pure-hydrogen WD simulation at T^ff = 5977 K and logg = 5.0 (top), and the 
stellar simulation at 6011 K, logg = 5.0, and [M/H] = —4 {bottom). 

from (3D) structures, hereafter (3D) spectra. For WD sim¬ 
ulations with log^ > 7.0, the normalized Balmer line profiles 
computed from (3D) structures are ne arly identical to the re¬ 
sults o f a full 3D spectral synthesis (iTremblav et al.l 1201 IL 
l2013cl) . However, we have found that the (3D) spectra ap¬ 
proximation becomes questionable at lower gravities, in line 
with main-sequence simulations where 3D inhomogeneities 
have a significant direct im pact on predicted Balmer line pro¬ 
files (iLudwig et al.l200^ . We review first the differences be¬ 
tween (3D) and ID spectra and then consider the additional 
effects from 3D spectral synthesis in Sect. 3.2. 

Figure |5] shows the comparison of predicted (3D) and ID 
Balmer line profiles for three simulations at logg = 6.0, which 
is an illustrative case for ELM WDs. The ID s pectra adopt the 
ML2/ a = 0.8 parameterization of the MET (ITremblav et al.l 
1201 Ol) and the same microphysics as the 3D calculations. The 
ELM results can be compar e d to th e case at logg = 8.0 in 
Eigure 16 of ITremblav et al.l (l2013al) . One major difference 
is that up to ten Balmer lines are typically fitted simultane¬ 
ously for ELM WDs, since the non-ideal effects and the sen- 
sitivity to surface gravi ty are shifted to higher Balmer lines 
(iGianninas et al.ll2014al) . Eor the simulation at Teff ~ 9000 K 
and logg = 6.0 in the middle panel of Eigure |5l it is only for 
HIO and higher lines in the series that (3D) profiles become 
shallower than ID profiles due to differences in the predicted 
non-ideal effects. The latter are very sensitive to the d ensity 
stratification in the atmosphere (ITremblav et al.ll2013cl) . The 



log Tr 


Figure 4. Temperature structure versus log tr for {3D) (red, solid) and ID 
ML2 /q = 0.8 (black, dashed) model atmospheres at ~7000 K and gravities 
from logg = 5.0 (bottom) to 9.0 (top) in steps of 0.5 dex. The temperature 
scale is correct for the logg = 5.0 model but other structures are shifted by 
2 kK relati ve to each other for clar ity. Properties of the logg > 7.0 models 
are given in ITremblav et alj <2013cr . 

shape of the lower Balmer lines is also significantly different 
between (3D) and ID spectra at intermediate Tj-ff, which is 
mostly responsible for Teff shifts. 

The (3D) and ID Balmer lines become identical at high 
Teff values (see the right panel of Eigure |5]|. Eor ELM WDs, 
this transition appears at relatively low temperatures, and 3D 
effects are largely negligible on the thermal structures for 
Teff > 10,000 K. Similarly to C/O-core WDs, convection be¬ 
comes very inefficient in this regime and the thermal structure 
is essentially fixed by the radiation field, even though convec¬ 
tive velocities are still large. Eor cool values, the left panel 
of Eigure|5]suggests that 3D effects are small, but Section 3.2 
demonstrates that the (3D) approximation is not always good 
in this regime. 

Eigure^ presents the (3D) atmospheric parameter correc¬ 
tions found by fitting the (3D) spectra with our standard grid 
of ID spectra. The corrections were derived simultaneously 
for the nine Balmer lines from H/3 to H12 in the same way 
we fit observations. The (3D) structures deviate significantly 
from their ID counterparts in the upper layers (tr < 10“^) due 
to the cooling effect of convective overshoot. This results in 
deep (3D) cores for Ha and H/3, as observed in the middle 
panel of Eigure |5] While we have no suggestion that the 3D 
simulations are inaccurate, there is no observational evidence 
for an adiabatic surface cooling in WDs and it appears prema¬ 
ture to account for these effects. As a consequence, the line 
cores were partially removed from t he fits, as well as t he entir e 
Ha line, as discussed in Sect. 3.1 of iTremblav et al.l (l2013cl) . 

The (31^ atmospheric parameter corrections for ELM WDs 
in Eigure |6] follow the trend observed for higher gravity ob¬ 
jects, despite the fact that more lines are included in the fitting 
procedure in the former case. The (3D) logg corrections have 
a similar strength for ELM and C/O-core WDs, although they 
are restricted to a narrower range of Teff at low surface gravity. 

3.2. 3D Spectral Synthesis 
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Figure 5. Compaiison of the blue wing of ten Balmer line profiles (Ha to H12) calculated from {3D) stnjctures (red, solid) and standard ID ML2/a = 0.8 
structures (black, dashed) for three models at logg = 6.0 from our sequence identified in Table[T] Teff values are shown on the different panels. All line profiles 
were normalized to a unity continuum at a fixed distance from the line center. The spectra were convolved with a Gaussian profile of 3 A resolution (FWHM) to 
represent typical observations. 
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Figure 6. (3D) atmospheric parameter corrections (solid black) found by 
fitting our grid of (3D) spectra with the reference grid of ID spectra relying 
on the ML2/a = 0.8 parameterization of the MLT. The ID = 3D reference 
parameters are on the intersection of the dotted lines, and 3D corrections 
are read by following the solid lines from the intersection. We utilized a 
resolution of 3 A and th e cores of the deeper lines were removed from the fits 
iTremblav et alJ[2013dl . The red dashed line represents the position of the 
maximum strength of H/3 for ID models. Tabulate d values are ava i lable in 
Table[3]and|4]for the range 5.0 < logg < 6.5 and in ITremblav et all 42013d) 
for the range 7.0 < logg < 9.0. The thin solid blue lines represent 1.5D 
atmospheric parameter corrections (see Section 3.2). 

We have employ ed the Linfor3D three-di mensional spec¬ 
tral synthesis code (iLudwig & Steffenll2008l) to compute 3D 
H/3 profiles for ELM white dwarfs. Unlike our previous ex¬ 
periments for log,g > 7.0, the 3D spectra and their (3D) coun¬ 
terparts from properly averaged (3D) structures differ by a 
few percent, with the largest effects at lower gravities. It is 
currently out of the scope of this work to proceed with the 
time-consuming computation of a full grid of 3D synthetic 


spectra. However, to further constrain the precision of the 
3D corrections, we have recomputed our grid of ELM spectr a 
with the so-called 1.5D approximation (Isteffen et alJll99-f) . 
Under this formalism, we assume that the physical conditions 
do not vary in the horizontal direction for each point on the 
top of the computational box, i.e. each column is a plane- 
parallel atmosphere. The 1.5D spectrum is then the average 
over all 150 x 150 columns, for which the emergent ID spec¬ 
tra can be easily calculated from our regular spectral synthesis 
code. The full 3D spectral synthesis has the effect of coupling 
nearby grid points, hence it is expected to lie somewhere be¬ 
tween the extreme cases of the 1.5D and (3D) approxima¬ 
tions. 

In Eigure|6l we compare the 1.5D (blue) and (3D) (black) 
atmospheric parameter corrections. We notice that while the 
(3D) approximation is relatively good at logg = 6.5, espe¬ 
cially in terms of logg corrections, it becomes questionable 
at lower gravities. We remind the reader that the 1.5D ap¬ 
proximation overestimates 3D effects, hence suggests an up¬ 
per limit to the 3D corrections. Since both 1.5D and (3D) 
corrections have the same signs, we suggest that applying the 
(3D) corrections provides a minimum 3D effect and a reason¬ 
able estimate of the atmospheric parameters of ELM WDs. 

The deficiency of the (3D) approximation for ELM WDs is 
likely explained by multiple factors. The temperature fluctua¬ 
tions are slightly larger at lower gravities although this could 
not alone explain the variations. On the other hand, we find 
that the line to continuum emergent intensity ratio has a sig¬ 
nificantly larger spatial variation for lower surface gravities. 
It implies that fluctuations are less likely to cancel in the 3D 
spectral synthesis. This could be related to the fact that broad¬ 
ening due to collisions with neutral atoms becomes increas¬ 
ingly important at the low Teff values where strong 3D effects 
are observed in ELM WDs. Einally, we have verified that 
the (3D) approximation is not the cause for the slight resid- 
ual bump in the 3D WD mass distribution at T^ff ^ 12,000 K 
(ITremblav et al.l 1201 3d iGenest-Beaulieu & BergerMl2014^ 
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Table 3 

ID ML2/a = 0.8 to 3D Tgff Corrections 


Table 4 

ID ML2/q; = 0.8 to 3D log^ Corrections 


Teff 

(K) 

logg = 5.0 
(cm s“^) 

logg = 5.5 
(cm s““) 

logg = 6.0 
(cm 

logg = 6.5 
(cm s“^) 

6000 

6 

12 

13 

24 

6500 

-33 

-19 

-11 

7 

7000 

-64 

-6 

5 

21 

7500 

-203 

-50 

-19 

16 

8000 

-460 

-202 

-105 

-66 

8500 

-555 

-425 

-197 

-152 

9000 

-288 

-564 

-338 

-154 

9500 

-99 

-243 

-557 

-151 

10000 

0 

-155 

-237 

-466 

10500 

0 

78 

-116 

-211 

11000 

0 

0 

136 

13 

11500 

0 

0 

0 

123 


Teff 

(K) 

logg = 5.0 
(cm s“^) 

logg = 5.5 
(cm s“^) 

logg = 6.0 
(cm 

logg = 6.5 
(cm s“^) 

6000 

0.010 

-0.019 

-0.017 

-0.018 

6500 

-0.023 

-0.037 

-0.040 

-0.024 

7000 

-0.128 

-0.117 

-0.122 

-0.134 

7500 

-0.239 

-0.225 

-0.215 

-0.239 

8000 

-0.218 

-0.316 

-0.315 

-0.350 

8500 

-0.064 

-0.181 

-0.302 

-0.340 

9000 

-0.035 

0.054 

-0.199 

-0.283 

9500 

0.001 

0.006 

0.056 

-0.150 

10000 

0.000 

0.044 

0.021 

0.004 

10500 

0.000 

0.037 

0.027 

0.002 

11000 

0.000 

0.000 

0.034 

0.007 

11500 

0.000 

0.000 

0.000 

0.024 


or the presence of deep cores for the lower Balmer lines com¬ 
pared to ID spectra. These features are also present when 
using 1.5D spectra. 

3.3. 3D Correction Functions 

We have derived functions to convert ID atmospheric pa¬ 
rameters of ELM WDs (logg < 7.0) to their 3D counter¬ 
parts. We have ensured that the transition is s mooth at log.? ^ 
7.0 wit h the functions already presented in iTremblav et alJ 
(I2013cll . which are more precise in the range 7.0 < logg < 
9.0. Eortran 77 functions are available in the online Appendix 
A. One can also rely on the data from Tables[3]and|4]to derive 
alternative functions or request (3D) spectra to the authors. 
Our independent variables are defined as 

gx = logg [cms“^]- 8.0 , ( 1 ) 


Tx = 


reff[K]-10,000 
1000 


( 2 ) 


and the fitting functions for the ML2/a = 0.8 parameterization 
of the MET are given below with the numerical coefficients 
identified in Table |5] 


ATe 


eff 


1000 


[ML2/a = 0.8]= a 13 -I- ^ai -I- < 34 exp 


(as + aq exp[fl8rx 


+ a<)gx\)Tx +aegx 


exp 


^(^ 2 + 


aioexp[aii7x-l-ai2gx])7x + a3.?x ), 


(3) 


Alogg [ML2/a = 0.8]= bu + (bx + biQe.x^[biiTx + bi 2 gx]) 
exp {[b 2 + b^ exp(/75 Tx + b(,gx}\Tx 

+ [bq, + bq Qxpib^Tx + ^9^x)]^x) ■ 

(4) 


Table 5 

Coefficients for 3D Correction Functions 


Coefficient 

(Teff) 

Value 

Coefficient 

(logg) 

Value 

ai 

-7.3801851E+00 

bi 

5.1729169E-02 

32 

6.0978389E-01 

b2 

-2.8514123E+00 

a3 

-8.5916775E-01 

b3 

-2.8616686E+00 

34 

7.2038240E+00 

b4 

3.2487645E+00 

35 

-1.1195542E-01 

bs 

-7.3122419E-02 

36 

-8.2204863E-03 

b6 

-3.5157643E-02 

37 

1.0550418E-01 

bv 

2.0585494E+00 

38 

2.2982103E-03 

bg 

2.9195288E-01 

39 

-2.3313912E-02 

b9 

-9.2196599E-02 

3l0 

-1.5513621E-01 

bio 

-3.2042870E-01 

3ll 

6.5429395E-01 

bii 

-6.3154984E-01 

312 

-1.1711317E+00 

bi2 

8.3885527E-01 

3|3 

2.8245633E-03 

bi3 

-1.0552187E-04 


4. ASTROPHYSICAL IMPLICATIONS 
4.1. The ELM WD Sample 

iGianninas et al.l (12014ah present a ID spectroscopic analy¬ 
sis of 61 low-mass WDs identified in the ELM Survey. All 
but five of these objects are in confirmed binary systems with 
periods ranging from 12 min to about a day. The majority of 
these targets are warmer than 12,000 K, where the 3D correc¬ 
tions are negligible. However, there are 23 convective WDs 
that are in the right Teff and logg range for 3D corrections to be 
significant. Table | 6 ] presents the updated parameters for these 
23 systems, including temperature, surface gravity, mass, ra¬ 
dius, cooling age, merger time, and their gravitational wave 
strain. The corrections are largest at around 8000 K for low- 
mass WDs. Hence, the parameters for relatively cool objects 
like J0345H-1748 (NETT 11748) have changed significantly. 
Since the 3D corrections lower the estimated logg, they also 
imply lower masses and gravitational wave strains for these 
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Table 6 

Convective ELM WD Parameters after 3D Corrections 


Object 

Teff 

(K) 

(cm s“^) 

M 

(M©) 

D 

(pc) 

R 

(O.OIR©) 

'^cool 

(Gyr) 

^2, min 
(Mo) 

^2,1=60° 

(M©) 

_max 
' merge 

(Gyr) 

logh 

JOl 12-1-1835 

9740 ± 140 

5.77 ± 0.05 

0.160 

662 ± 98 

8.63 ± 1.07 

1.8 

0.62 ± 0.03 

0.85 ± 0.04 

2.68 

-22.40 

J0152-I-0749 

10800 ± 190 

5.93 ± 0.05 

0.168 

980 ± 131 

7.36 ±0.83 

1.4 

0.57 ± 0.03 

0.78 ± 0.04 


-22.80 

J0345-I-1748 

8560 ± 120 

6.51 ±0.04 

0.181 

175 ± 34 

3.91 ±0.41 

3.9 

0.76 ± 0.03 


7.18 

-22.16 

J0745-I-1949 

8230 ± 130 

5.91 ± 0.07 

0.156 

300 ± 56 

7.28 ± 1.07 

3.6 

0.10 ±0.01 

0.12 ±0.02 

5.80 

-22.65 

JOS 18-1-3536 

10190 ± 190 

5.87 ± 0.09 

0.164 

3388 ± 640 

7.76 ± 1.26 

1.7 

0.25 ± 0.04 

0.33 ± 0.06 

9.30 

-23.47 

J0822-I-2753 

8800 ± 130 

6.36 ± 0.05 

0.175 

496 ± 72 

4.57 ±0.51 

4.2 

0.76 ± 0.08 

1.05 ±0.11 

8.18 

-22.31 

J0849-I-0445 

10020 ± 150 

6.27 ± 0.06 

0.176 

1101 ± 158 

5.08 ± 0.61 

1.9 

0.65 ± 0.04 

0.89 ± 0.05 

0.45 

-22.38 

J0900-I-0234 

8330 ± 130 

5.96 ± 0.07 

0.155 

618 ± 113 

6.83 ± 0.99 

3.8 





J1005-1-3550 

9760 ± 140 

6.00 ± 0.05 

0.166 

1231 ±392 

6.72 ± 0.79 

2.1 

0.19 ±0.02 

0.24 ± 0.02 

10.63 

-23.12 

J1112-I-1117 

9240 ± 140 

6.17 ±0.06 

0.169 

266 ± 42 

5.57 ± 0.70 

2.7 

0.14 ±0.01 

0.17 ±0.02 

12.73 

-22.56 

J1443-I-1509 

8970 ± 130 

6.44 ± 0.06 

0.181 

557 ± 85 

4.24 ±0.51 

3.5 

0.84 ±0.12 

1.17 ±0.17 

3.82 

-22.25 

J1512-I-2615 

11250 ± 180 

6.93 ± 0.06 

0.250 

766 ± 96 

2.85 ±0.31 

0.4 

0.31 ±0.04 

0.39 ± 0.05 


-22.94 

J1518-I-0658 

9650 ± 140 

6.68 ± 0.05 

0.197 

309 ± 41 

3.38 ± 0.35 

0.9 

0.58 ± 0.03 

0.78 ± 0.04 


-22.42 

J1538-I-0252 

10030 ± 150 

5.98 ± 0.06 

0.167 

1154 ± 180 

6.90 ± 0.92 

1.9 

0.76 ± 0.06 

1.06 ±0.09 


-22.86 

J1557-1-2823 

12470 ± 200 

7.76 ± 0.05 

0.448 

191 ± 18 

1.45 ±0.11 

0.4 

0.41 ± 0.03 

0.51 ±0.04 


-21.92 

J1614-I-1912 

8700 ± 170 

6.32 ±0.13 

0.172 

207 ± 54 

4.73 ± 1.00 

4.1 





J1630-1-2712 

10200 ± 170 

6.03 ± 0.08 

0.169 

2161 ± 374 

6.56 ± 0.99 

1.8 

0.52 ± 0.04 

0.70 ± 0.06 


-23.13 

J1741-I-6526 

10410 ± 170 

6.02 ± 0.06 

0.170 

998 ± 150 

6.71 ±0.87 

1.7 

1.11 ±0.05 

1.57 ±0.06 

0.16 

-22.11 

J1840-1-6423 

9120 ± 140 

6.34 ± 0.05 

0.177 

740 ± 109 

4.70 ± 0.54 

3.3 

0.65 ± 0.03 

0.88 ± 0.04 

4.78 

-22.47 

J2103-0027 

9900 ± 150 

5.79 ± 0.05 

0.161 

1253 ± 177 

8.47 ± 1.01 

1.8 

0.71 ±0.04 

0.98 ± 0.05 

5.73 

-22.72 

J2119-0018 

9690 ± 150 

5.72 ± 0.08 

0.158 

2835 ± 507 

9.09 ± 1.38 

1.7 

0.74 ± 0.04 

1.03 ±0.05 

0.58 

-22.82 

J2228-I-3623 

7890 ± 120 

5.78 ± 0.08 

0.175 

420 ± 79 

8.96 ± 1.34 

7.7 





J2236-I-2232 

11310 ± 170 

6.54 ± 0.04 

0.182 

355 ± 43 

3.81 ±0.40 

1.1 

0.39 ± 0.03 

0.51 ±0.04 


-22.79 


Notes. J0345+1748 is NLTT 11748 and J2236+2232 is LP 400-22. 


systems. _ 

ID spectroscopic analysis of iGianninas et al.l (1201 4al) find 
Teff = 8680 ± 120 K and log^ = 6.83 ± 0.04 for NLTT 1 174i|. 
The ID to 3D corrections from our models revise these to 
Teff = 8560± 120 K and logg = 6.51 ±0.04, a change of 0.3 
dex in surface gravity. Using the lAlthaus et alJ (1201 3h evolu¬ 
tionary models, the revised mass, distance, and radius are M = 
0.18±0.02MQ,£/=175±34pc,andR = 0.0391±0.0041RQ, 
respectively. Our distance estimate is in excellent agreement 
with the para l lax me asured at the USNO, tt = 5.6 ±0.9 mas. 
iKaplan et alj (1201 4|) use high-quality photometry of NLTT 
11748 to constrain the same parameters from the observed 
eclipses. Depending on the thickness of the surface hydro¬ 
gen layer of the secondary (more massive) WD, they derive 
Ml =0.136-0. 162M0 andRi =0.0423-0. 0433Rq. Our mass 
and radius measurements, including the 3D correction factors, 
are consistent with their results within the errors. 

4.2. The Mass-Radius Relation 

There are several tidally distorted and/or eclipsing ELM 
WD systems that provide model-independent constraints on 
the radii. NLTT 11748 is discussed above; it is an eclipsing 
system where the spectroscopically inferred radius from the 
ID models were significantly smaller than the eclipse mod - 
elling suggested (iKaplan et al.ll20141 : IGianninas et ai.ll2014a^ . 
Studying the radius co nstraints on the ELM WD sample, 
IGianninas et al.l (1201 4al) find that the agreement between the 
spectroscopically inferred radius and the model-independent 
values is quite good for Teff > 10,000 K. However, there 

’ Tlie earlier analysis by IKilic et all j2010bll found logg = 6.54 ±0.05. 
However, this was ba sed on older WD models wit hout the improved Stark 
broadening profiles of ITremblav & BergeronI 1200911 . 


are four relatively cool objects, LP 400-22, NLTT 11748, 
J0745, and J0751, where the radius is underestimated or the 
mass/logg is overestimated. Ligure [T] shows the comparison 
of the updated radius measuremen ts, including the 3D cor - 
rections, for the same 10 stars from iGianninas et^ (l2014ah . 
Table |2] provides the references for the different radius mea¬ 
surements. The agreement between the spectroscopic and 
photometric radii is now excellent over the temperature range 
8,000-17,000 K. There are only two d eviant points, LP 400- 
22 and J0745, and we refer the reader to lKilic et al.l ( 120131) and 
IGianninas et al.l (l2014bh . respectively, for further discussion 
of these objects. Ignoring LP 400-22 and J0745, the excellent 
agreement between the radius measurements from these two 
independent methods suggests that the 3D model correction 
factors presented in this paper have the right amplitude to ac¬ 
count for the differences seen in the earlier analysis relying on 
ID models. Linally, we note that this independent verification 
of 3D results is u nique to ELM WDs and fairly different to the 
one presented in iTremblav et al.l (l2013d) for single C/O-core 
WDs with trigonometric parallax measurements. 

4.3. ELM WDs with Milli-Second Pulsar Companions 

There are several WD companions to milli-second pulsars 
with optical spectroscopy available in the literature. These 
WDs provide an independent method to estimate the mass 
of the neutron star through constraints on the mass of the 
WD and the semi-amplitude of the radial velocity variations. 
Hence, accurate mass (and logg) measurements for the ELM 
WDs in these systems are crucial. Many of these WDs have 
temperatures in the range 8000-10,000 K, where 3D effects 
are significant. 

PSR J0218±4232, J1012±5307, and J1909-3744 have 
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Table 7 

ELM WDs with Photometric Radius Measurements 


Object 

Rff 

(K) 

Rspec 

(O.OlRo) 

^phot 

(O.OIR©) 

Ref. 

J0056-0611 

12230 ± 180 

5.65 ±0.61 

5.6 ±0.6 

1 

JO106-1000 

16970 ± 260 

6.42 ± 0.68 

6.3 ±0.8 

1 

J0112-H835 

9740 ± 140 

8.64 ± 1.07 

8.8 ±0.9 

1 

J0345-H748 

8560 ± 120 

3.91 ±0.41 

4.23 ±0.10 

2 

J065U-2844 

16340 ± 260 

3.25 ±0.31 

4.0 ± 0.2 

1 

J0745-H949 

8320 ± 130 

7.28 ± 1.07 

17.6f|» 

3 

J0751-0141 

15740 ± 250 

13.15 ± 1.04 

13.8^:? 

1 

J174U-6526 

10410 ± 170 

6.71 ±0.87 

7.6 ±0.6 

1 

J2119-0018 

9690 ± 150 

9.09 ± 1.38 

10.3 ± 1.6 

1 

J2236-b2232 

11310 ± 170 

3.81 ± 0.40 

>9.9 

4 


References. ( DIHermes et alj j2014i) : f2) IKaplan et aljf20T4l) : (3) this work; 
(4) |Kilic_ekJI (2013D. Spe ctroscopic parameters are taken from this work 
and IGianninas et a TI I2014an. _ 

Notes. See also iGianmn^ et al.l 42014bl) regarding the i?phot measurement 
uncertainties for J0745+1949. 
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Figure 7. Differences in radii as measured from our 3D spectroscopic anal¬ 
ysis and photometric analyses using eclipses or ellipsoidal valuations. The 
error bars represent the en'ors of the two independent radius measurements 
added in quadrature. T he errors for J0745+1949 reported in Figure 11 of 
IGianninas et alJ i2014ah are incorrect. 


companions with Teff ~ 8000-85 00 K and logjg < 7 . based 
on ID spectroscop i c analyses dvanKeAwijket^ 1199^ 
ICallanan et al.l 119981: iBassa et aP l2003t LTacobv et al.ll2005l) . 
The corrections in log^ are on the order of 0.3 dex at t hese 
temperatures. Eor example, Ivan Kerkwiik et ahl (Il996h de¬ 
rive Teff = 8550 ± 25 K and logg = 6.75 ± 0.07 for the ELM 
WD companion to PSR J1012-I-5307. Taking the 3D correc¬ 
tions into account brings down these values to Teff = 8440 K 
and logg = 6.43. PSR 1911-5958A has a WD companion 
with T pff = 10090 ± 150 K and logg = 6.44 ± 0.20 (iBassa et al.l 
1 20061) . The amplitude of the 3D corrections is lower at 
these temperatures, bringing down the temperature and sur¬ 
face gravity to 9780 K and logg = 6.38, respectively. 

PSR J1738-I-0333 is an excellent system where the mass. 


radius, and the surface gravity of the WD companion can 
be estimated from bo th optical spectro s copy/p hotometry and 
radio observations. lAntoniadis et al.l (l201^ obtain Teff = 
9130±150 K, logg = 6.55 ±0.10, M = O.181!jj;[’g^M0, and 
R = O.O37 ;|^qq3 R0 from optical spectroscopy. Interestingly, 
they also derive logg = 6.45 ±0.07 (from the orbital period 
decay rate of the binary, mass ratio, parallax, and photom¬ 
etry), M = 0.182 ±O.O16M0 (from the mass ratio and orbital 
decay rate of the binary), and R = 0.042 ±0.004/^0 (from pho¬ 
tometry). The latter measurements are independent of the ID 
spectroscopic analysis and imply a lower surface gravity and 
a larger radius for this ELM WD. These provide additional 
evidence for the “high logg problem” in ID spectroscopic 
analysis of ELM WDs, though (admittedly) the errors in sur¬ 
face gravity and radius are relatively large. Based on the 3D 
corrections factors that we estimate, we revise the spectro¬ 
scopic parameters of this companion to Teff = 8910 ± 150 K 
and logg = 6.30 ±0.10, which are consistent (within the er¬ 
rors) with logg = 6.45 ±0.07 estimated using the independent 
method. We not e that this WD co mpanion also happens to be 
a pulsating WD (iKilic et al.ll2015h . 

4.4. The ZZ Ceti Instability Strip 

Eigure |8] shows the instability strip for the pulsating DAV 
WDs based on ID and 3D spectroscopic analyses. Here we 
only show the par ameters for the DAVs derived from a uni¬ 
form analysis by IGianninas et al.l (1201 IL l2014al) . We also 
show the five ELM pulsators currently known (excluding the 
companion to PSR J1738-I-0333 since the ID model spectro¬ 
scopic analysis of that star was not performed using the same 
models), and a large number of photometrically constant stars 
that help define the boundaries of the instability strip. The 
solid lines show the boundaries of the instability strip from 
Gianninas et al. (2014). 

There is a clear trend of increasing logg with decreasing 
temperature in the ID spectroscopic parameters of these stars. 
This trend disappears in the 3D version of the plot, where the 
3D corrections lower the implied surface gravity. The blue 
edge of the instability strip is well defined in the ID analysis. 
The 3D corrections move the stars systematically to slightly 
cooler temperatures and lower surface gravities. The dashed 
lines show the revised boundaries based on the 3D corrected 
parameters. Eor the blue edge, the boundary is simply shifted 
to the right by 100 K. The blue edge of the instability strip is 
still well defined for both normal mass and ELM WDs in the 
3D analysis. Eor the red edge, as in the ID case, we use GD 
518 and SDSS J222859.93-1-362359.6 as the two anchor points 
at either extremity (the highest and lowest surface gravity) to 
define the 3D red edge, which now has a higher slope. The 
equations of the revised blue and red edges are defined by 
Equations |3 and |6] respectively. Given the small number of 
ELM pulsators currently known, it is too early to definitively 
constrain the red edge of the instability strip for those stars. 
The new ZZ Ceti instability strip from the 3D analysis does 
not reveal any new ELM pulsator candidates in the published 
list of ELM WDs, but the revised boundaries should be used 
in the future to look for new pulsating DA WDs. 

(logg)biue = 5.96923 X 10-4(reff)biue + 0.52431 (5) 

(logg)red = 8.06630 X 10-4(7;ff),ed-0.53039 (6) 

5. CONCLUSIONS 
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Figure 8. The ZZ Ceti instability strip based on the ID (left panel) and 3D (right panel) spectroscopic analysis. White diamonds represent the pulsators, whereas 
black dots show the ph otometrically constant WDs. The solid blue and red lines represent the empirical boundaries of the ZZ Ceti instability as determined by 
IGianninas et alj I2014af) . The dashed lines denote tentative boundaries matching the location of both instability strips (for normal mass and ELM pulsators) based 
on the 3D corrected parameters. The green error bars denote the remaining ELM WDs that have not yet been investigated for photometric variability. 


We provide mean 3D model spectra for WDs with 5 < 
log^ < 7, extending the previously published 3D model grid 
to very low-mass WDs. Our 3D models now cover the entire 
temperature and surface gravity range for the observed pop¬ 
ulation of convective DA WDs. We list ID to 3D correction 
factors, which can be as large as 0.35 dex in surface grav¬ 
ity. These corrections largely resolve the discrepancies seen 
in the logg and radius estimates for relatively cool ELM WDs 
that are observed in double WD or WD -i- milli-second pulsar 
systems. We also update the boundaries of the ZZ Ceti insta¬ 
bility strip for WDs, including the ELM pulsators as well as 
the normal (logg ~ 8) WDs. 
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APPENDIX 

CORRECTION FUNCTIONS FOR FORTRAN 77 

ID ML2/a = 0.8 to 3D T^ff corrections 


1 function ML18_to_3D_dTeff_ELM(Teff, logg) 

2 real*8 ML18_to_3D_dTeff_ELM, a(13) , Teff, logg 

3 real*8 TX,GX,Shift 

4 c 

5 c ; IN: Teff, effective temperature (K) 

6 c ; logg, surface gravity (g in cm/s2) 

7 c 

8 c ; OUT: Teff correction (K) ID ML2/alpha = 0.8 to 3D 

9 c 
10 

11 A(l) = 

12 A(2) = 

13 A(3) = 

14 A(4) = 

15 A(5) = 

16 A(6) = 

17 A(7) = 

18 A(8) = 

19 A(9) = 

20 A (10 ) 

21 A(ll) 

22 A(12) 

23 A(13) 

24 

25 TX=(Teff-10000.0)/lOOO.00 

26 GX=(logg-8.0) 

27 Shift=A(13)+(A(1)+A(4)*exp((A(5)+A(7)*exp(A(8) 

28 * *TX+A(9)*GX))*TX+A(6)*GX))*exp((A(2) 

29 * +A(10)*exp(A(ll)*TX+A(12)*GX))*TX+A(3)*GX) 

30 

31 ML18_to_3D_dTeff_ELM=Shift*1000.0 

32 return 

33 end 


-7.3801851E-f00 
6.0978389E-01 
-8.5916775E-01 
7.2038240E-f00 
-1.1195542E-01 
-8.2204863E-03 
1.0550418E-01 
2.2982103E-03 
-2.3313912E-02 
=-l.5513621E-01 
= 6.5429395E-01 
=-l. 1711317E-f00 
= 2.8245633E-03 


ID ML2/a = 0.8 to 3D log,g corrections 

1 function ML18_to_3D_dlogg_ELM(Teff,logg) 

2 real*8 ML18_to_3D_dlogg_ELM, B(13) , Teff, logg 

3 real*8 TX,GX,Shift 

4 c 

5 c ; IN: Teff, effective temperature (K) 

6 c ; logg, surface gravity (g in cm/s2) 

7 c ; 

8 c ; OUT: log g correction (g in cm/s2) ID ML2/alpha = 0.8 to 3D 

9 c ; 

10 
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11 

B(l) = 

5.1729169E-02 

12 

B(2) = 

-2.8514123E+00 

13 

B(3) = 

-2.8616686E+00 

14 

B (4) = 

3.2487645E+00 

15 

B(5) = 

-7.3122419E-02 

16 

B(6) = 

-3.5157643E-02 

17 

B(7) = 

2.0585494E+00 

18 

B(8) = 

2.9195288E-01 

19 

B(9) = 

-9.2196599E-02 

20 

B(10) 

=-3.2042870E-01 

21 

B(ll) 

=-6.3154984E-01 

22 

B(12) 

= 8.3885527E-01 

23 

B(13) 

=-l.0552187E-04 

24 



25 

TX=(Teff-10000.0)/lOOO . 0 

26 

GX=(logg-8.0) 

27 

Shift 

=B (13) + (B (1) +B (10) *exp (B (11) *TX+B (12) * 

28 

* 

*exp((B(2)+B(4)*exp(B(5)*TX+B(6)*GX)) 

29 

■k 

*TX+(B(3)+B(7)*exp(B(8)*TX+B(9)*GX))* 

30 



31 

ML18_ 

to_3D_dlogg_ELM=Shift 

32 

return 

33 

end 




